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Abstract. Quantum gravity may lead to an energy de- 
pendence in the speed of light. The high energy radiation 
from gamma-ray pulsars can be used to place limits on 
such effects. We find that emission from the Crab pulsar 
at energies above 2 GeV trails that at 70-100 MeV by 
no more than 0.35 ms (95% confidence) and place a lower 
bound on the energy scale of quantum gravitational effects 
on the speed of light of 1.8 x 10 15 GeV. This bound might 
be improved by two orders of magnitude by observation of 
pulsations from the Crab at higher energies, 50-100 GeV, 
in the near future. 
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1. Introduction 

Quantum gravity may cause modification of the disper- 
sion relation for photons at high energies. It has recently 
been suggested that certain quantum gravity models may 
lead to a first order correction to the dispersion relation 
which can be parameterized as At = LAE/cEqg, where 
At is the magnitude of the travel time difference between 
two photons whose energies differ by AE and that have 
traveled a distance L, and Eqq is the energy sca le of 
the dispersion effects ( Amelino-Camclia ct al. 1998). To 
probe dispersion effects at high energy scales, accurate rel- 
ative timing of nearly simultaneously produced photons 
of different energies which have traveled long distances 
is required. Use of sub-millisecond time structure of the 
keV ph oton flux of gamma-ray bursts at cosmologica l dis- 
tances ( Amclino-Camelia et al. 1998| ; Schacfcr 1998 ) and 
use of several minute time structure in TeV flares from 



AGN ( Biller et al. 1998 ) have been suggested. Here, we 
show that sub-millisecond timing of GeV emission from 
gamma-ray pulsars may also place useful constraints on 
the dispersion relation for photons at high energies. 

Below, we use existing gamma-ray data to determine 
the accuracy with which high energy pulsar emission can 
be timed and to place bounds on the energy scale for quan- 
tum gravity corrections to the speed of light. We then 



discuss how this limit might be improved by pulsar obser- 
vations at higher energies in the near future. 

2. Gamma-ray pulsations from the Crab 

We chose to analyze the energy dependence of pulse ar- 
rival times from the Crab pulsar as it has the largest ratio 
of distance to pulse period of the bright gamma-ray pul- 
sars, thus maximizing the constraints which can be placed. 
Also, the pulses from the Crab are well aligned in time 
from radio waves, through optical and x-ray emission, to 
gamma-rays. Thus, it is likely that the photons of different 
energies are produced nearly simultaneously. 

We used data from the Energetic Gamma-Ray Ex- 



periment Telescope (EGRET) (Thompson ct al. 1993) of 
the Compton Gamma-Ray Observatory (CGRO). We ex- 
tracted gamma-ray photon event lists from the CGRO 
public archive for observations pointed within 40° of the 
Crab and then, using the program pulsar (version 3.2, 
available from the CGRO Science Support Center) se- 
lected events lying within the energy dependent 68% point 
spread function of EGRET and calculated the phase of 
each photon relative to the radio ephemeris of the Crab 



(Arzoumanian et al. 1992). The pulse period of the Crab 



changed from 33.39 ms to 33.49 ms over the course of 
these observations. The radio timing must be corrected 
for the variable dispersion along the line of sight to the 
Crab. T he accurac y of this correction is estimated to be 
0.2 ms ( Nice 1998| ), consistent with previous estimates of 
the accuracy of the dispersion correction for the Crab 
( iGullahorn et al. 1977| ). 

Pulse phase histograms for several energy bands are 
shown in Fig. |l|. The main pulse peak, near phase 0.0, is 
the most appropriate feature for timing. The main peak 
i s similar acr oss the energy range from 70 MeV to 2 GeV 
( Ficrro 1995 ). The peak width is about 0.05 in phase, and 
appears somewhat narrower at high energies. There is no 
obvious shift of the peak centroid with energy. 

To study the energy dependence of the speed of light, 
we measured the main peak pulse arrival time in each en- 
ergy band. We did this in two ways. First, we calculated 
the average arrival time for photons in the main peak. We 
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Fig. 1. Crab pulsar phase histograms for various 7-ray 
energy bands. Zero phase is set by the radio ephemeris. 

found the average time for each energy band using photons 
with phases between -0.0464 and 0.0336, an interval cen- 
tered on the mean arrival time for all photons used in this 
analysis. Second, we parameterized the pulse arrival times 
by fitting a Lorentzian to the pulse profile, within the same 
phase range specified above, for each energy band. Before 
fitting, a constant rate equal to the average rate between 
phases -0.4 and -0.2 was subtracted. The resultant was 
then fit with a Lorentzian using a gradient-expansion al- 
gorithm to compute a non- linear least squares fit. The fits 
were all acceptable with x 2 m the range 2.9 to 7.6 for 5 
degrees of freedom. 

Fig. H shows the pulse arrival times calculated via both 
methods. The errors in Fig. |^ correspond to Ax 2 — 1 (68% 
confidence). The energy of each point is the median photon 
energy for each energy band. For the highest energy band, 
the median energy is substantially lower than the average, 
2.9 GeV versus 5.0 GeV. The zero pulse phase is set by 
the radio ephemeris. The pulse arrival time for all photons 
used in this analysis is shown as a dashed line and differs 
by 0.21 ms from the radio zero phase. T his is withi n the 
error in the radio dispersion correction ( Nice 1998|) . We 
note that errors in the radio zero phase can broaden the 
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Fig. 2. Pulse arrival time versus energy for the Crab. The 
diamonds indicate the average arrival time for photons 
within the main pulse for each energy band. The squares 
indicate the centroid of a Lorentzian fit to the pulse profile 
for each energy band. The energies plotted are the median 
energy for each band; the diamonds are shifted slightly in 
energy for clarity. The dashed line is the centroid of a 
Lorentzian fit to the pulse profile for all energies above 
70 MeV. 



gamma-ray peak, but will not induce an energy dependent 
shift in the gamma-ray pulse arrival time. The accuracy 
of the pulse arrival time determination for the Lorentzian 
fit is 0.07 ms (Ax 2 = 3.84 or 95% confidence for a sin- 
gle parameter of interest) in the 100-200 MeV band and 
0.21 ms (95% confidence) in the highest energy band. The 
accuracy in the highest energy band is limited mainly by 
statistics. 

It is apparent from the figure that there is no statisti- 
cally significant variation in pulse arrival time with energy. 
To place an upper bound on any energy dependence in the 
speed of light, we compare the arrival time for photons 
with energies above 2 GeV (median energy 2.93 GeV) to 
that for the 70-100 MeV band (median energy 82.8 MeV). 
The 95% confidence upper limit on the difference of the 
arrival time s is 0.35 ms. A dopting a distance to the Crab 
of 2.2 kpc ( Zombeck 1990 ), this leads to a lower limit on 
the energy scale of quantum gravity effects on the speed 
of light of E QG > 1.8 x 10 15 GeV (95% confidence). This 
limit lies below the range of interest, but within an order 
of magnitude of some predictions in the context of string 
theory (|Witten 1996|). 



3. Discussion 

Other effects which could also produce an energy depen- 
dent delay in photon arrival times include energy depen- 
dent dispersion due to the strong gravitational field near 
the neutron star, purely electromagnetic dispersion, an en- 
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ergy dependence in the emission location, or an intrinsic 
energy dependence in the emission time. The effect of any 
energy dependent dispersion due to the strong gravita- 
tional field near the neutron star is likely to be small be- 
cause, even if emitted from the neutron star surface, pho- 
tons traverse the region of high gravitational fields within 
about 0.1 ms. Allowing a fractional change in the speed 
of light equal to the dimensionless field strength at the 
neutron star surface, GM/Rc 2 w 0.2, where M « 1.4 M Q 
is the neutron star mass and R ss 10 km is the neutron 
star radius, the difference in arrival times would be only 
0.02 ms. The actual energy dependent change in the speed 
of light is likely to be much smaller than 0.2. Any signif- 
icant purely electromagnetic dispersion at MeV energies 
and above can be excluded based on the dispersions mea- 
sured at lower energies. 

An energy dependence in the photon emission loca- 
tion or intrinsic emission time could produce a significant 
energy dependent time delay. While the possibility that 
precise tuning of the emission locations or times for vari- 
ous energy photons could cancel an energy dependent dis- 
persion arising from quantum gravity effects, we consider 
such a coincidence unlikely, although not excluded, and 
interpret our lack of detection of any energy dependence 
in arrival times as constraining both the energy dependent 
dispersion and the emission location and time. In this case, 
the average emission location, projected along our line of 
sight, for photons at energies in the 70-100 MeV band 
must lie within 110 km of that for photons above 2 GeV, 
within 50 km of that for 0.5-1.0 GeV photons, and within 
150 km of that for radio photons. 

It is nnrniirap-ino; that flip analysis gtinws that it is 



der correction to the dispersion relation for electromag- 
netic waves, then measurement of the pulse arrival time 
of the Crab at 50 GeV with an accuracy of 0.1 ms could 
be used to place a lower bound on Eqq > 1.1 x 10 17 GeV. 
This is within the range, 10 16 — 10 18 GeV, for the energy 
scale for quantum gravity effects preferred in string theory 



(Witten 1996) 



If future measurements do reveal an energy depen- 
dence in pulsar photon arrival times, then it will be dif- 
ficult to distinguish an energy dependent dispersion from 
an intrinsic energy dependence in the emission location 
or emission time. This problem is common to all of the 
suggested astronomical tests of quantum gravity effects. 
Convincing proof for quantum gravity effects will likely 
require detection of energy dependent time delays in at 
least two different classes of objects, preferably at vastly 
difference distances, i.e. pulsars versus AGN or gamma- 
ray bursts, with all of the detections compatible with the 
same value of Eqq . 
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